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We present a detailed asymptotic analysis of correlation functions for the two com- 
ponent spanning tree on the two-dimensional lattice when one component contains 
three paths connecting vicinities of two fixed lattice sites at large distance s apart. 
We extend the known result for correlations on the plane to the case of the upper 
half-plane with closed and open boundary conditions. We found asymptotics of cor- 
relations for distance r from the boundary to one of the fixed lattice sites for the 
cases r ^ s ^ 1 and s 7' 1. 
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INTRODUCTION 



In recent years the logarithmic conformal field theories ( 



m 



CFT) [a] and their relation to 
lattice models of statistical physics like dense polymers js], y , , the sandpile model g], 3, S Q] , 
dimer models 10| and percolation j^, have been the subject of active research. Among 
all these models, the Abelian sandpile model (ASM) [l^ is one of the most interesting 
and fruitful, because correlation functions containing logarithmic corrections can be found 
explicitly using combinatorial methods. In this way, the full correspondence between the 
lattice model and the logarithmic conformal field theory can be transparently tested. A 



lot of successful checks have been made in 



including various calculations for 



correlations in the bulk and on the boundaries, the determination of boundary changing 
fields, the insertion of isolated dissipation and evaluation of some finite-size effects. One of 
the most popular questions and checks was about the origin of the logarithmic corrections 



2 



to the correlations. It is known just two cause of that: the insertion of the dissipation at 
isolated sites 3], that follows from the logarithmic behavior of the inverse Laplacian at a 
large distance, and non-locality of height variables for h > 2 in the ASM joj. The last 
case is more important and difficult, because sandpile configurations with height variables 
h > 2 are mapped onto an infinite set of non-local configurations of spanning trees. This 
non-locality arises due to the presence of a specific three-leg subgraph, so-called 6-graph 
[igI . The 0-graph is a sub configuration of the spanning tree, consisting of three paths, that 
connect the vicinity of vertex jo with that of vertex (FigH]). The paths with additional 
branches attached to them form one component, which is surrounded by another component 
of the spanning tree. A generalization of 0-graph is an odd "fc-leg" subgraph, which has 
been considered by E.V.Ivashkevich and C.-K.Hu in [l^ for the infinite square lattice. They 
obtained the asymptotic dependence P(r) ~ Inr/r 2 , for k = 1,3,5,... and concluded 
that it is the presence of the second component is responsible for the logarithmic correction 
to the correlation function. Indeed, H.Saleur and B.Duplantier considered correlations of 
"k-leg operators" by mapping of the two-dimensional percolation problem on a Coulomb gas 
and found, that the asymptotics of correlation functions for one component spanning tree 
has a pure power-law decay r^^^~2^ [l^ . 

Later on, G. Piroux and P. Ruelle calculated the height probabilities for /i > 2 in the ASM 
on the upper-half plane with closed and open boundaries 9|. They enumerated the spanning 
tree configurations with 0-graph, having one fixed site at distance r apart the boundary 
and another site running over the whole upper half-plane. The summation over positions 
of the running site leads to cumbersome estimations of integrals, so that it is difficult to 
follow details of correlations between different parts of G-graph explicitly. Calculations of 
two point correlations Pih — PiPh in the ASM on the plane for h > 2 also lead to the 
same difficulties. The G-graph arising for h > 2 consists of three paths connecting one 
fixed site with height h > 2 with another site at the distance s, running over the whole 
plane except the site with the height variable "one". It was shown in 15|, that evaluation 
of the logarithmic corrections to the two-point correlations does not need the summation 
by s over the whole plane. Instead, it is enough to take into account only those G-graphs 
configurations when the running site of the G-graph is situated in a vicinity of the height 
variable "one". The latter approach, being much simpler, is not so transparent, and an 
additional analysis of correlations of different parts of the G-graph is desirable. In this 
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work, we find the asymptotic behavior of three-leg correlations for the case of the upper-half 
plane. We test validity of the method described in for the upper half-plane, examining 
the order of expansion, where we can obtain a disagreement. 

II. THE MODEL 

We consider the labeled graph Q = (V, E) with vertex set V and set of bonds E. The 
vertices are sites of the square lattice and an additional point which is the root 'V: V = 
{sx,y, {x,y) G Z^, |x| < M, \y\ < N} U {-k}. The bonds of E connect only neighboring sites. 
Vertices ixi,yi, jx2,y2 ^ ^ neighbors, if {x2 — Xi)^ + (2/2 — Ui)^ = 1- Also all boundary 
vertices {s^^±Ar,x G Z} and {s±M,y,y ^ Z} are neighbors of the root -k. The graph Q 
represents a finite square lattice of 2N + 1 x 2A/+ 1. In thermodynamical limit, N,M ^ 00, 
the lattice Sx,y covers the whole two dimensional plane. We consider also the upper-half plane 
with closed and open boundary conditions at the lattice sites Vq = {sx,i,x G Z, |x| < M}. 
The corresponding graphs are Qd = and Qop = (V+,-Eop), where V+ = {sx,y,x G 

Z,y eN, \x\ < M,y < N}U {i<}, with the sets of bonds E^ = E D {{i,]), i G I4, j G V+} 
and Eop = E^i U {(z,*), i G Vq}. We construct the desired spanning tree configurations on 
the above graphs by using the arrow representation, see e.g. Accordingly, we attach 

to each vertex i G an arrow directed along one of bonds G E incident to it. 

Each arrow defines a directed bond {i — > i') and each configuration of arrows ^ on ^ defines 
a spanning directed graph (digraph) Qdiri-^) with set of bonds Edir{A) = i G 

V"\{^}, i' G V, G E} depending on A. Similarly, the arrow configurations on Q^i 

and Qop define a spanning digraph Gci,dir{A), Gop,dir{A) with corresponding sets of bonds. 
Note that no arrow is attached to the root so that it has out-degree zero. A sequence 
of directed bonds (zi, 12), (^2, "^3), (^3, "^4), • • • , {im-i, im) is called the path of length m from 
the site ii to the site im- This path forms a loop if im = ii- Spanning tree is a spanning 
digraph without any loops. Our aim is to construct a two-component spanning tree, with 
one component containing three paths connecting neighboring sites jo, jo — x,jo + y with 
sites to, to — y,to + x, where jo and to have coordinates (0, r) and {k,l + r) respectively, 
and X = (1,0), y = (0,1) are unit vectors (FigH]). The relevant configurations will be 
investigated with aid of the determinant expansion of the Laplace matrix. This technique 

nnn 

is described in details in [9|, [12, 116|. Let the vertices of the set V, be labeled in arbitrary 




FIG. 1: Three- leg subgraph: the simplified version of G-graph. 



order from 1 to = = (2M + 1)(2A^ + 1). Then Laplacian A of size n x n has the 

elements: 

Zi iii = j, 

— 1 if i,j are nearest neighbors, (1) 
otherwise. 

where Zi is the degree of vertex i G ^\{*}. The determinant of Laplace matrix is equal to 
the number of spanning trees on Q with the root -k. Laplace matrices Aop, Ad for the upper 
half-plane have size n = = (2M + 1)A^ and are defined by the same way as A, but 

for graphs Qop, Qci respectively. The determinant of A is a sum over all permutations a of 
the set {1,2,..., n}: 



det A = ^ sgn(o-) Ai,<^(i) A2,a(2) • • • A„,„(„ 



(2) 



o-GS„ 



where Sn is the symmetric group, sgn = ±1 is the signature of permutation a. In general, 
each permutation a E Sn can be factorized into a composition of disjoint cyclic permutations, 
say, cr = Ci o C2 o . . . Cfc. This representation partitions the set of vertices ^\{*} into non- 
empty disjoint subsets which are orbits O = {vi i,Vi^2 ■ ■ ■ j'^iJ,} C V of the corresponding 
cycles Ci, i = 1, . . . , k, at that U^^^Oi = V"\{^} and Yli=i — where k is the length of 
cycle Cj. The orbits consisting of just one element, if any, constitute the set Sfp{a) of fixed 
points of the permutation: Sfp{a) = {v = a{y),v E ^\{^}}- A cycle q of length \ci\ = k > 2 
is called a proper cycle. The proper cycles on Q are of even length only, hence, the number 
of proper cycles p defines the signature of the permutation a, that is sgn(o") = (—1)^. Thus 
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Eq. 



can be written as follows: 



n [n/2] 

detA = Y[z, + J2i 

i=l p=l 



n 



a=c^ o...oCu i=l 



{vi),vi 11 



(3) 



where is the /c-fold composition of the cyclic permutation q of even length li, Vi E Oi(cr), 
so that c\~^{yi) ^ c^{vi) and c['(i;j) = Vi. The term YYi=i^i equals to the number of all 
spanning digraphs Qdir{A), having the root -k. Each of others terms on the right-hand side 
of Eq.(l3]) having a non-zero set of fixed points Sfp 7^ up to a sign equals to nje5^p(o-) 
because all non-diagonal elements equal to —1. That product represents the number of 
distinct spanning digraphs which have in common the specified cycles Ci, . . . ,Cp, and differ 
in the oriented edges outgoing from vertices j G Sfp{a). These oriented edges may form 
cycles on their own which do not enter into the list Ci, . . . , Cp. The proper cycles formed by 
the oriented bonds incident to fixed points of a given permutation a = ci o C2 o . . . o Cp should 
enter into enlarged list of cycles Ci o C2 o . . . o Cp o . . . o Cp' , p' > p, corresponding to another 
permutation a'. The expansion ([3]) can be interpreted in form of the inclusion- exclusion 
principle [l6|. Let Ci, C2, . . . , be the list of all possible proper cycles on Q. We define Ai, 
2 = 1, 2, . . . , m as the set of all spanning digraphs on Q containing the particular cycle q 
and Aq is the set of all spanning digraphs Qdiri-^)- Let Ast be the set of spanning trees on 
Q. Then we can write of Eq.(l3]) in the form: 



det A = \A 



ST I 



\An 



l<i<j<m 



n A 



(4) 



1=1 



where |y4| is cardinality of the set A. E g .(HI) is the Kirhhoff theorem for the number of 
spanning tree subgraphs of a given graph jl7 |. 



III. THREE-LEG CORRELATIONS 



Now we modify the Laplace matrix changing three non-diagonal elements: 
Zi if z = j, 

— 1 if i,j are nearest neighbors, 

-e if (z, j) G i3 = {(jo,to), {jo-x,to-y), {jo + y,to + x)} 
otherwise. 



(5) 
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FIG. 2: Component I consists of tree paths and branches attached to them, component II is another 
spanning tree surrounding component I. 

In the same determinant expansion as Eq.([3]), only the terms containing the product 
Ajy Ajy_£^to-yAjo+j/,io+£ ^ survive in the hmit hm^^oo det A'/e^. Permutations a = 

Ci o C2 o . . . o Cp corresponding to these terms, contain cycles Cp^ o . . . o Cp^ with directed bonds 
B. Since sites from B form angles on the lattice, topologically we can draw only one or three 
cycle(s) containing these bonds. Thus, expression lime^oo det A' /e^ equals to the number of 
configurations with following features: (i) each configuration is a two component spanning 
graph on the plane; (ii) one component consists of three paths connecting the vicinity of 
site jo with the vicinity of site to and branches of the spanning tree attached to these paths; 
(iii) another component is the spanning tree having the root ^, and surrounding the first 
component. The quotient of such configurations and all one-component spanning trees is: 

det A' , detfJ + ^G) , 
hm - — -— — - = - hm = det A 6 

e^oo (-£:)3det A e^oc ^ ^ ^ ^ 

where 5 = A' — A and G = A^^ = G[(ni,mi); (n2,m2)] is the Green function, which is 
defined in thermodynamical limit M, TV — > oo on the plane as: 

G[{ni,mi);{n2,m2)]^^^^^ = G(ni - na, mi - ms) = 

= Go,o+ // ^ , (7) 

J J-T, oTT^ 2 — COS a — cos p 

In Appendix we give more details about this function, including its asymptotics on a long 
distance for arbitrary direction of the vector between (ni,mi) and {n2,'m2). For the case of 
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the upper half-plane, presence of the boundary changes the Green function: 

Gd = A~i^ [{rii, mi); (n2, 1712)] = G{ni - 712, mi - 7712) + G{ni - n2, rui + m2 - 1) (8) 
Gop = A~p [(rii, mi); (^2, m2)] = G{ni - n2, mi - m2) - G{ni - n2, mi + m2) (9) 

Matrices in Eq. ([6]) are of size n x n (or n+ x n+ for the upper half-plane), but, matrix 5 
has only three non-zero elements: 



v 



Jo Jo -X jo + y 

-e 

-£ 

-e 



to 

to + X 



(10) 



so, we obtain due to ([6]) the matrices of size 3x3 only: 





Gk,l 


Gk+i,i 


Gk,i~i 


\ 




Gk,l~l 


Gk+i,i-i 


Gk,l-~2 






\^Gk+i,i 


Gk+2,l 


Gk+i,i^ 


V 



(11) 





Gk,i - 


^ Gk,l+2r~l 


Gk+i,i - 


- Gk+lJ+2r~l 


Gk,l^l 


+ GkJ+2r 


Ad = 


Gk,i~i 


+ Gk,l+2r~2 


Gk+i,i~i 


+ Gk+l,l+2r~ 


2 Gk,l~2 H 


~ Gk,l+2r~l 




yGk+l,l - 


^ Gk+l,l+2r- 


1 Gk+2,i - 


- Gk+2,l+2r~l 


Gk+l,l-l 


+ Gk+l,l+2r j 




Gk,i - 


~ Gk,l+2r 


Gk+i,i — 


Gk+l,l+2r 


Gk,i-i — 


r \ 

tJfc,Z+2r+l 


Aop 


Gk,i-i - 


~ Gk,l+2r-l 


Gk+i,i-i — 


Gk+l,l+2r-l 


Gk,i-2 - 


~ Gk,l+2r 




\Gk+i^i - 


" Gk+l,l+2r 


Gk+2,i — 


Gk+2,l+2r 


Gk+i,i~i — 


Gk+l,l+2r+l j 



(12) 



(13) 



For further analysis we find it convenient to split determinant expansions into four groups: 

detA,Kfc,/,r) = F(=''°)(A;,/) -F(''^)(fc,/,z) +F(1'2)(A:,/,z) -F(°'=^)(fc,z) (14) 
detAop(A;,/,r) = F^^^''\k,l) ^ F^''-^\k,l,z) ^ F^^-^\k,l,z) ^ F^''^^\k,z) (15) 

where 2; = Z + 2r for open boundary and 2; = Z + 2r — 1 for closed one. 
Functions F^^'"^ are: 



F^'''\kJ)=detApiane{k,l) 



(16) 
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F^^'^\k, /, z) = —Gk,zGl_^i i^i + Gk,zGk,i-2Gk+2,i + G k^z-iG k+i,iG k+i,i~i — Gk,z-iGk,i~iGk+2,i 
— '^Gk+i,zGk+i,iGk,i~2 + "^Gk+i^zGk^i-iGk+i^i-i — Gk+i,z~iGk,iGk+i,i-i + 
+ Gk+2,zGk,iGk,i~2 — Gk+2,zG\i_i + Gk+i,z~iGk+i,iGk,i-i — 



Gk+i,z+iGk,iGk+i^i~i — Gk,zGk,iGk+2,i + Gk+i,z+iGk+i,iGk,i~i 



k,z+l^k,l~l^k+2,l — ^k,z^k+l,l + ^ k,z+l^ k+l,l^ k+1, 



(17) 



F^^''^\k,l, z) = —Gk,zGk+2,zGk,i~2 + Gk+i,z~iGk,zGk+i^i-i — Gk,z-iGk+i,zGk+i,i^i + 

+ Gk,z-iGk+2,zGk,i^i + G\j^i.,Gk,i-2 — G k+i,z-iG k+i,zG k,i~i — Gk,zGk+2,zGk,i + 
+ Gk,zGk+i,z+iGk+i,i-i + G\.^Gk+2,i — Gk,z-iGk+i,z+iGk+i,i + 2Gk+i,zGk,zGk+i,i + 
+ Gk,z-iGk,z+iGk+2,i — Gk+i,zGk,z+iGk+i,i-i + Gk ,z+lGk+2,zGk,l-l + 
+ Gk+i,z-iGk+i,z+iGk,i — G k+i,zG k+i,z+iG k,i-i — Gk,zGk+i,z-iGk+i,i (18) 



F^^'^\k, z) = —Gk,zGk+i,z-iGk+i,z+i + G\^^Gk+2,z + Gk,z-iGk+i,zGkJri,z+i 



(19) 




F"'-=''(k,z) F<'-2'(k,l,z) 

FIG. 3: Angle-angle system near boundary. 



The meaning of functions F'^^'^^ is shown in FiglHl The index /i = 1, 2, 3 shows numbers 
of "short" hnks between the vicinity of site jo and the vicinity of site to and index v = 1,2,3 
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shows number of "long" links between the vicinity of the mirror image Jq of site Jq and the 
vicinity of to- 
ll sing the asymptotic formula for the Green function (Appendix, ( 1501) ) we can analyze 
behavior of functions F^^''^^ for s = |s| ^ 1 and \s + d\ ^ 1, where s = {k,l), d = (0, rf); 
d = 2r for the case of UHP with the open boundary and d = 2r — 1 for the closed one: 



F(3'0)(A;,/) 

F^^'^\k,l,l + d) 
F^^'^\k,l,l + d) 
F^'''^\k,l + d) 



1 1 + In s + Co 



+ 



1 2cos^/3 — 1 1 cos a cos/? 1 ln|s + (i|+co 



+ 



1 2 cos^ a — 1 ^ 1 cos a cos /5 ^ 1 In s + cq 

87r3 |^+ J|4 



+ 



s^\s + d\^ 47r3 |s + rf|3s 
1 1 + ln|s + J| + Co 



+ . . . 



in-" 



\s + d\' 



+ ... 



(20) 

(21) 

(22) 
(23) 



where a and (3 are angles between horizontal axis and vectors s, d + s correspondingly. The 
constant is cq = — 27rG'o,o + 7 + | In 2. 




(-1) 




FIG. 4: For large d S> s 1 two configurations are canceled due to asymptotical similarity of 
paths. 

The function F^^''^\k, 1,1 + d) has the leading term l/d^ for d ^ s ^ 1, and its contribu- 
tion tends to zero as (i — > oo. This function describes configurations with two "long" links 
and one "short" link. We can see from FigJH that such configurations define two types of 
spanning tree subgraphs with opposite sign (this sign changed because the number of loops 
containing "long" links is changed by one). Paths on the square lattice in these two types 
of subgraphs are topologically almost identical, and differ only in the vicinity of the site J'q. 
The difference disappears, when d oo. 
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Now we consider expansions (l20!) - (l23l) of (fT^ . (fT5|) in two ways: first we expand these 
expressions by r assuming r ^ s ^ 1; second we expand them by s assuming s ^ r ^ 1. 
In both cases we take only a = or a = 7r/2 for simphcity. As a result, we get for r s ^ 1, 
a = 0: 



det Ad(s,r) 



1 + ln(2r) + In s + 2co 



O 



o 



1 



3 2 + ln(2r) + In s + 2co ^ § + ln(2r) + In s + 2co ^ ^ 



2s4 



+ 



(24) 



and for a 



det Kiis.r 



1 + ln(2r) + Ins + 2co 



1 



IQ-n'^r \ s 



+ 



1 1 + ln(2r) + In s + 2co f| + ln(2r) + In s + 2co 



24s2 6s3 

For s ^ r ^ 1 and a = we obtain: 

3 + 2co + 2 In s 



2s4 



O 3 (25) 



+ 



det Ad(s,?^ 
2 



r — r I + 
10 



(co + In s)(6r^ - 12r^ + + 5r) + — (3r^ - 6r^ + + 2r) + O — (26) 



and for a 



2 • 



det Ki{s,r) 



+ 



-^(l + 2co + 21ns 

1 



2(co + lns)(5r^-4r) + (3r^-2r) +0 



Asymptotics of det Aop(s,r) for open boundary conditions for r s ^ 1, a = is 



(27) 



det Aop(s,r) 



1 - ln(2r) +lns 



+ 4 



+ 



1 / 5 ln(2r)-lns f + 3 ln(2r) - 3 In s 



^ + 



IGvr^r^ \ 4s^ 



0'^ 



oil 



(28 



For a = I we have 



det Aop(s,r) 



1 - ln(2r) + Ins [I 



Svr^s^ 



1 



3r \ 2s3 s^ 



+ 



3 i-ln(2r)+lns if + 3 ln(2r) - 3 In s /I 



IGyr^r^ \ 4s2 



(29) 
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In the range of s ^ r ^ 1 we have for a = 0: 



det Aop(s, r) 



and for a 



det Aop(s,r) 



12 



4(r2 + 3r3 + 2r^) 



74r^ - 394r6 + 416rS - 96ri° ^ / 1 

+ — 



14 



37r3s9 



15 



+ 



10 



(30) 



(31) 



IV. DISCUSSION 



100 





-0.0010 



FIG. 5: det Aop(s,r) and AetKd{s,r) for s = {k,l), r = 50, -50 < A: < 50, -50 < / < 50. 



As it was noticed above, the 6-graph is the key object in calculations of heights variables 
hi > 2 in the ASM. The previous analysis 0, Q, [s, 9] shows that ASM belongs to a c = —2 
minimal model of logarithmic conformal field theory. Height variable h = 1 is associated 
with a primary field (f){z,^) with conformal weights (1, 1) and heights /i = 2, 3, 4 behave 
like its logarithmic-partner iplz,!) with scaling dimension 2. Let Py^^^{r) be the probability 
of height h = 1,2,3,4 at distance r apart from the boundary of the upper half-plane. In 
fact this is a two point correlation function, and by mixing operators, its dependence on r 
enables to obtain a structured constant for the operator product expansion (OPE) and to 
conjecture the logarithmic behavior of two-point correlations Whi,hj{f) = Phi,hj{f') — PhiPhj 
of height variables /ij > 1 on the plane at site i and height variables hj > 2 at site j at 
distance r apart from the site j. Correlations P^^^(r) have been obtained combinatorially, 
by mapping ASM onto the spanning trees model j^. 

Due to similarity between the 0-graph and the three-legs correlations, we may compare 



calculations in with the present ones. The "head" of 0-graph located at site i with hi = 2 
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corresponds to one of the fixed points of the three-legs correlations, the running point of 
0-graph corresponds to another point at distance s apart. The dependence of det A(s,r) 
on s for fixed r is shown in Figj5] for open and closed boundary conditions. We see that 
the function has a strong peak at {k,l) = (0,0). An essential part of computational work 



m [9| is summation over all positions of the running point. Instead, we may try to use the 
expansions (l20l) - (!23l) for integration in the vicinity of the peak using the fact that det A(s, r) 
decays as Ins/s^ with s. 

In the case of two-point correlations Wi,2('"); this method leads to drastic simplification of 
calculations due to rapid convergence of the integrals over the vicinity s < sq of the peak for 
So << r [l^. In the case of boundary correlations P2^^{r), the leading term of asymptotics 
by r in (125|) and ( 129|) is In r/r^ and its integration by s in the vicinity of the peak is not 
sufficient for obtaining a coefficient at In r/r^. Indeed, the first terms of expansions fl25|) and 
(!29|) contain Ins/s^ which gives also In r/r^ upon summation over the half-plane: 

oo oo , I -,1 , 

EE^i^^-^^. (32) 



But expansions (l25il and ( l29l) are not valid for s ~ r and therefore the method 15|] fails in 
the case of calculations of P2^^{r). 

Despite the failure with the description of the 0-graph near the boundary, the expansion 
(!24|) - (l3T]) are still useful for a LCFT treatment. The previous attempts to describe the 
logarithmic correlations both in the lattice and field theories were undertaken solely for the 
ASM. In this case, the logarithmic partner of the primary field is associated with the height 
variables hi > 1 having a non-local representation in the spanning tree model. Moreover, 
the non-local representation is the infinite sum over positions of the running point of the 
G-graph. But two branching points of the 6-graph are in turn some correlating objects of 
the spanning tree which can be considered in the framework of the LCFT independently of 
the ASM problem. Thus the collection of asymptotics (124|) - (13T]) for the three- leg correlations 
near closed and open boundaries, together with the bulk asymptotics (1201) should be found 
within the LCFT provided that one finds a proper identification for the three-leg branching 
points. 
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V. APPENDIX 

We consider here the Green function Gp ^ for two-dimensional square lattice and derive its 
asymptotic expansion for large distances = ^Jp^ + ^ 1 j using the methods proposed 
in and in 18]. We denote an angle between the vector r = (p, q) and the horizontal axis 
by V9, so that 

p=rcosv9, (33) 
q=r sin ip. (34) 

Since the Green function contains singular part Gq.Oj we consider a function ^ = Gp g— Gq.Oj 
which has an integral representation 

1 f ^ipa+iql3 _ ^ 

9p,q = ^ o n^a d/3. (35) 

ovr^ J 2 — COS a — cos p 

and obeys the symmetry relations, 

9p,g — 9-p,q — 9p-q — 9q,P) (36) 

so we can put p > g > without loss of generality. After the integration over a and 
symmetrization by (3 we come to the expression 

1 /-TT I 2 — cos 13 - ^J {2 - cos - l) cosqp - 1 
27r Jo v^(2 - cos/3)2 - 1 

Consider the Taylor expansion of the function 

^ (2 - cos /? - V(2-cos/3)2-l) ^ (38) 

for small positive up to order 14 and denote it Q(/3): 

0(B) = _ ^ , (112p=^-493p)/j9 

^^^^ 12 96 288 40320 1152 1161216 

421p2/5io (I27741p-46200p3)/5ii (140/ - 3887^^) 

^ 1935360 ^ 1277337600 ^ 69672960 ^ 

(69432p3 _ 152461p) p^'^ (629861^^ - 43120/) p^'^ 

^ 6131220480 ^ 42918543360 ' ^ 

The function 

R^{p) = (^2 - COS /5 - V(2-cos/5)2 - 1)^ - e-P^Q{p) (40) 
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vanishes at /3 = 0, and has a unique maximum over (3 for all p. The location of the maximum, 

can be found as a series in ^ if we construct the series expansion of the derivative '^^^^^ 

and recursively equate coefficients to 0. Calculations give 

15 ^ 14503977 

28p3 + 17248p5 ^ ■ ■ ■ ^ 



and 



It is easy to show that 



^ 120135498046875 ^, 



3"e-'''cos(,/3)^^~^, n>0 (43) 



It means that we can replace the function in the integral (l37j) by its Taylor expansion in (5 
up to 14-th order and get an expression for the Green function with accuracy 0(p~^°). It is 
convenient to express the Green function as a sum of three integrals 



where 



9p,q = /l + /2 + /3, (44) 



i,=Lr(t:i^^\u,B (45) 



, e-P^QM cos(g/3) e'^'P cos{q/3) \ dp ^ ^ f 1 , 
''-I ' V(2-cos^P-l P )2^ + ^l;^)- 







The expressions give 

Ji = (Inr + 7 + In vr) + exponentially small terms (48) 

27r 

1 / 31n2\ , , 

Finally we obtain 

_lnr + 7+^ cos(4(^) 1 /3cos(4(^) 5cos(8(^) 
^^'^ ~ 2^^ ^ 247rr2 ^ ^ \ 8(h^ ^ 96^? 

1 /51cos(8(^) 35cos(12(^)\ 
^ ^ V 224^^^ ^ 144^^ ) ^ 

1 / 217cos(8yp) , 45cos(12y.) , 1925 cos(16y.) ^^ 1 \ 

+ 1^ 6407r + 16^^ + 768^^ ^ ^ [p^ J ^^"^ 
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where the expansion 



1 



1 p 43/33 949/35 21727/3^ 501451/3^ 



+ 



COS/3)2 - 1 



13 12 1440 120960 9676800 766402560 



+ 



8112267073/3^ 



277899049/3^3 
4782351974400 



+ o ip'') 



(51) 



41845579776000 



is used for the calculation of fjlTl) . 
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